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I.  INTRODUCTION 


The  use  of  separate  codes  for  computing  inviscid  flow  and  turbulent 
boundary  layer  development  over  yawed,  spinning  and  non-spinning  bodies  of 
revolution  has  yielded  some  very  good  solutions  for  cone  and  ogive-cylinder 
shapes1 2 3 4 5.  However,  the  authors  of  Reference  1  have  found  that  application  of 
these  techniques  ,to  bodies  with  boattailed  afterbodies  has  not  yielded 
satisfactory  results  even  at  small  angles  of  attack  (a  <  4°). 

Several  recent  publications  have  reported  supersonic  flow  field  computa¬ 
tions  using  Parabolized  Navier-Stokes  (PNS)  techniques.  These  publications 
have  reported  very  good  results  for  cone  models  for  laminar  and  turbulent 
viscous  flow  *  ,  and  for  cone  and  ogi ve-cyl inder  models  for  laminar  viscous 
flow  .  The  PNS  method  appears  to  offer  an  attractive  technique  for 
computing  flow  oyer  bodies  with  discontinuities  in  surface  curvature  (such  as 
occurs  at  the  junction  between  the  cylinder  and  the  boattail)  since  the 
inviscid  flow  and  viscous  layer  are  computed  simultaneously.  Further,  the  PNS 
method  permits  adequate  flow  field  resolution  to  be  achieved  with  very 
reasonable  computer  costs.  This  report  documents  detailed  comparisons  of  PNS 
computational  results  to  experimental  measurements  of  turbulent  boundary  layer 
profile  characteristics  of  a  spinning  ogi ve-cyl inder- boattail  body  at  Mach  = 
3.  In  addition,  comparisons  are  made  between  the  PNS  computations,  boundary 
layer-invi scid  computations,  and  experimental  measurements  of  Magnus  forces 
for  cone,  ogive-cylinder,  and  ogi  ve-cyl  i  nder-boattai  1  models  for  2  <  M  <  4. 
The  PNS  code  is  that  reported  by  Schiff  and  Steger6  .  The  boundary  layer- 
invi  scid  code  is  that  reported  by  Sturek ,  et  al 1  . 


1.  Stuvek ,  W.B.,  Dwyer,  H.A.,  Kaysev,  L.D.,  Nietubicz,  C.J.,  Reklis,  R.P., 

and  Opalka,  K.O.,  " Computations  of  Magnus  Effects  for  a  Yawed,  Spinning 
Body  of  Revolution,  "  AIM  Journal,  Vol.  16,  No.  7,  July  1978,  pp.  687-692. 

2.  Lin,'  T.  C.,  and  Rubin,  S.  G.,  "Viscous  Flow  Over  a  Cone  at  Moderate 

Incidence:  I  Hypersonic  Tip  Region",  Computers  and  Fluids,  Vol.  1,  1973 
pp.  37-57.  - 

3.  Lubard,  S.  C .,  and  Helliwell,  W.  S.,  "Calculation  of  the  Flow  on  a  Cone  at 

High  Angle  of  Attack",  AIAA  Journal,  Vol.  12,  July  1974,  pp.  965-974. 

4.  Rakich,  </.  V.,  Vigneron,  Y.  C . ,  and  Agarwal,  R.,  "Computation  of 
Supersonic  Viscous  Flows  Over  Ogive-Cylinders  at  Angle  of  Attack",  AIAA 
Paper  No.  79-0131,  17th  Aerospace  Sciences  Meeting,  January  1979. 

5.  Agarwal,  R.,  and  Rakich,  J.  V.,  "Computation  of  Hypersonic  Laminar  Viscous 
Flow  Past  Spinning  Sharp  and  Blunt  Cones  at  High  Angle  of  Attack",  AIAA 
Paper  No.  78-65,  AIAA  16th  Aerospace  Sciences  Meeting,  January  1978. 

6.  Schiff,  L.  B.,  and  Steger,  J.L.,  "Numerical  Simulation  of  Steady 

Supersonic  Viscous  Flow",  AIAA  Journal,  Vol.  18,  No.  12,  December  1980 
pp.  1421-1430.  - 
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II.  OVERVIEW  OF  NUMERICAL  SCHEME 


A*  Governing  Equations  and  Numerical  Scheme 

A  body-conforming  n,  ?,  coordinate  system  (Figure  1)  is  used  which 
maps  the  body  surface  and  outer  boundary  of  the  flow  region  in  physical  space 
onto  coordinate  surfaces  of  the  computational  space.  This  transformation 
simplifies,  the  application  of  surface  boundary  conditions  and  permits  the 
approximation  of  neglecting  streamwise  and  circumferential  viscous  terms  in 
hi  gh-Re.ynolds-number  flow  (see  Ref.  6).  The  resulting  steady  thin- layer  PNS 
equations  can  be  written  in  strong  conservation-law  form  in  terms  of  nondimen- 
sional  variables  as 


3F  a  a  ^ 

_ls  +  JF  ,  _9G  _  1  as 
95  3c  "  ^  a c 


(1) 


where 

£  =  £(x)  is  the  axial  (marching)  coordinate 

=  ^(XjyjZ)  is  the  circumferential  coordinate 

5  =  s(x,y,z)  is  the  normal  coordinate 

The  inviscid  flux  vectors  in  Eq.  (1)  are 


pU  “ 

P  V  “ 

puU+?xP 

,  F  =  J"1 

PuV+nvp 

A 

pvU 

pvV+n  p 

pwU 

J 

PwV+nzp 

(e+p$)U 

1 

(e+p)V 

L  J 

pW 

PuW+?xp 

pvW+CyP 

pwW+czp 
(e+p) W 


(2) 


with  contravari ant  velocity  components 

U  "  «xu 

v  =  nxu  +  HyV  +  nzw  L 
w  =  Cxu  +  eyv  +  Czw 


(3) 


8 


The  internal  energy  of  the  gas  e^  is  defined  in  terms  of  the  conservative 
variables  as 

e.j  =  (e/p)  -  0.5(u2  +  v2  +  w2)  (4) 

while  the  equation  of  state  for  a  perfect  gas  with  ratio  of  specific  heats  y 
i  s 

P/p  -  (Y  -  l)ei  =  a2/y  (5) 


Variations  of  body  geometry  are  included  in  Eq.  (1)  through  the  presence  of 
the  metric  terms  Sx,  nx,  ny,  etc.,  which  appear  in  the  flux  vectors.  The 

thin-layer  viscous  terms,  valid  for  hi gh-Reynolds-number  flow,  are 


S 


y(?H+?zK  +  (»‘/3)(CxUe+Cyv5+C2w5)5x 
u(Vcy+?z)v4  +  (w/3)UxVSVczwc)cy 

w(?x+?y+?z)w4  +  ^3H  WWV^Z 
{ ( ?x+ey+cz) t (u/2) (u2+v2+w2) ^ 

+  KPr'^Y-irV2)^  +  (u/3)(cxu 
+  5yv  +  ?zw)(cxu5  +  4yV?  +  4zw?)} 


(6) 


Equation  (1)  is  parabolic  with  respect  to  £,  and  can  thus  be  marched  downstream 
in  the  %  direction  from  an  initial  data  plane  (subject  to  appropriate  body  and 
free-stream  boundary  conditions)  under  those  conditions  where  the  local  flow 

A 

is  supersonic.  By  evaluating  the  pressure,  ps,  in  the  E$  flux  vector  using 

the  subsonic  layer  approximation,  Eq.  (1)  can  be  kept  stable  for  marching  for 
subsonic  points  as  well.  When  ps  is  set  equal  to  the  local  pressure  for 

supersonic  points,  and  is  evaluated  from  3p$/a4  =  0  (Figure  2)  for  points 

within  the  subsonic  viscous  layer  adjacent  to  a  wall,  Eq.  (1)  can  be  marched 
stably  for  all  flows  where  U  >  0;  that  is,  for  flows  without  streamwise 
reversal . 

The  numerical  algorithm  used  to  march  Eq.  (1)  downstream  is  an 
approximately-factored,  fully  implicit,  finite-difference  scheme.  The 
algorithm  can  be  written  in  so-called  delta  form  as 
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[A$J  +  (1  -  a)A?(6TiBj)](Asj)"1 

X  [K  J  +  (1  -  a)AS(6.CJ  -  Re'1  6  MJ)  |AqJ' 
s  4  5  J 


=  _(A$J  -  A^_1)qJ  +  a(E$J  -  E^"1) 

-  (1  -  oc)A?{6n[n^+1(E/J)J  +  ny+1(F/J) J  +  nj+1(G/J)J' 

+  6J?x+1(E/J)J  +  ^+1(F/J)j  +  4+1(G/J)J] 

-  Re~ll^}  -  [(€x/J)j+1Epj  -  (CX/J)JEJ-1]  +  Vqi 


(7) 


A  -1 

where  q  =  J  (p,  pu,  pv,  pw,  e). 

The  6' s  represent  second  order  central  difference  operators  while  A 
represents  a  conventional  forward  difference.  The  Jacobian  matrices  A,  B.  and 

3E  3F  dC  9  9 

c  are  defined  as— ,  and^j,  respectively.  The  coefficient  matrix  M  is 

obtained  from  the  Taylor  series  linearization  of  the  viscous  vector  S.  The 
algorithm  shown  in  Eq.  (7)  is  second  order  accurate  in  C  for  a  =  1/3,  and  is 
first  order  accurate  in  5  for  a  =  0.  The  fourth  order  dissipation  term  is 
represented  by  V  which  is  added  to  damp  high-frequency  oscillations. 

The  alqorithm  is  conservative  and  of  second-order  accuracy  in  the 
marching  direction.  A  two-layer,  Cebeci-type  eddy  viscosity  model7  is 
included  for  the  computation  of  turbulent  flows.  The  algorithm  has  been 
applied  to  compute  a  variety  of  laminar  and  turbulent  flows  and  the  results 
have  been  in  excellent  agreement  with  those  obtained  from  more  costly  time- 
dependent  computations.  Full  details  of  the  Parabolized  Navier-Stokes 
assumption,  and  of  the  derivation  of  the  algorithm  are  included  in  Reference 
6. 


B.  Conical  Initial  Solutions 


In  general  the  initial  data  plane  for  the  marching  method  must  be 
supplied  from  an  auxiliary  computation.  However,  when  treating  the  flow  over 
conical  or  pointed  bodies,  the  marching  code  can  be  used  to  generate  its  own 
initial  data  plane.  As  outlined  in  Ref.  6,  for  inviscid  conical  flows  a 


7.  Baldwin ,  B.  S and  Lomax,  H.,  "Thin  Layer  A pp roximati on  and  Algebraic 
Model  for  Separated  Turbulent  Flows",  AIAA  Paper  No.  78-257,  16th 
Aerospace  Sciences  Meeting,  January  1978. 
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conical  grid  is  selected  and  the  flow  variables  are  initially  set  to  free- 
stream  values.  The  solution  is  marched  downstream  from  an  initial  station 
and,  after  each  step,  the  solution  is  scaled  to  place  it  back  at  the  original 
station.  When  no  change  in  the  flow  variables  occur  with  further  marching, 
the  variables  are  constant  along  rays,  and  a  conical  solution  has  been 
generated.  If  the  flow  variables  within  the  viscous  layer  can  also  be  assumed 
to  be  locally  constant  along  rays,  the  same  procedure  can  be  used  to  generate 
viscous  conical  solutions. 

For  the  ogi ve-cyl inder-boattai  1  computations,  the  tip  of  the  ogive  was 
replaced  with  a  cone  tangent  to  the  ogive  at  x  =  0.267  caliber  (1  caliber  =  1 
model  diameter).  Conical  solutions  were  generated  at  that  station  and  used  as 
starting  data  for  the  marching  code.  A  small  error  is  made  in  generating  the 
starting  solution  for  cases  with  surface  spin  due  to  the  change  in 
ci rcumferential  velocity  with  longitudinal  position.  However,  this  error  is 
small  since  the  stepsize  used  in  the  conical  initial  solution  is  less  than  1% 
cf  the  distance  from  the  tip  to  the  initial  plane.  In  any  event,  this  error 
is  quickly  corrected  as  the  initial  solution  is  marched  over  the  body. 


C.  Adaptive  Grid 

An  adaptive  grid  capability  was  developed  for  the  PNS  code  in  order  to 
maintain  adequate  resolution  of  the  viscous  layer  as  the  solution  develops 
over  the  full  length  of  the  spinning  model.  This  capability  was  also  found  to 
be  of  importance  in  maintaining  computational  stability,  particularly  for  the 
Mach  =  4  computations.  The  strategy  used  here  was  to  check  for  the  value  of 

y+  ^y+  =  pwUTy//uw  ’  Ut  =  at  the  flrst  9rid  node  above  the  model  sur¬ 

face  and  to  adjust  the  grid  stretching  parameters  to  maintain  this  value  of  y+ 

within  the  desired  range,  5  <  y+  <  10.  This  check  was  made  only  at  the  wind 
and  lee  sides  of  the  model.  The  stretching  parameter  was  varied  linearly 
between  the  extremes  determined  at  the  wind  and  lee  sides  for  grid  nodes  at 

circumferential  stations  off  the  pitch  plane.  If  the  value  of  y+  was  found  to 
be  outside  the  specified  criteria,  the  stretching  parameter  was  adjusted  by 
0.5%  for  the  next  computational  step.  Although  the  adaptive  grid  technique 
works  well  in  general,  it  was  found  that  permitting  too  great  a  change  in  the 
grid  configuration  from  step  to  step  resulted  in  the  introduction  of  errors 
into  the  computation. 


III.  RESULTS 

A*  Model  Geometry  and  Experimental  Measurements 

The  dimensions  of  the  ogi ve-cyl inder-boattai 1  model  used  for  the  detailed 
f  ow  field  studies  are  shown  in  Figure  3.  The  model  is  6  calibers  long  with  a 

1-caliber,  7  boattail,  and  closely  resembles  a  modern  low-drag  artillery 
projecti le. 
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A  number  of  wind-tunnel  experiments  have  been  conducted  for  this  model 
geometry  in  order  to  obtain  data  for  comparison  to  numerical  computations. 

The  data  acquired  include  measurements  of  wall  static  pressure8,  turbulent 
boundary- layer  velocity  profiles9’10,  surface  skin  friction9  ,  and  flow 

visualization.  The  test  conditions  were  M  =  3  with  a  tunnel  total  pressure  of 
0.298  MPa  and  tunnel  total  temperature  of  308°K.  These  conditions  produced  a 
free-stream  Reynolds  number  of  7.3  x  106  based  on  the  model  length.  The 
boundary  layer  was  tripped  near  the  tip  of  the  model  to  produce  a  reliable 

turbulent  flow.  Additionally,  aerodynamic  force* 11’ ^  measurements  were 
available  for  10°  cone,  ogive-cylinder  and  ogi ve-cyl i nder-boattai 1  shapes  at 
M  -  2,  3  and  4.  These  data  include  Magnus  and  pitch  plane  aerodynamic  force 
and  moment  coefficients.  All  tests  were  performed  using  SSWT  Number  One  at 
the  U.S.  Army  Ballistic  Research  Laboratory.  This  facility,  which  is  no 

longer  in  operation,  was  a  continuous  flow  tunnel  with  a  flexible  plate  noz¬ 
zle.  The  test  section  size  was  330  x  380mm  (13  x  15  in.). 


B.  Comparison  Between  Computation  and  Experiment 

Computations  were  performed  for  a  body  having  the  same  geometric  shape  as 
the  experimental  model,  and  for  flow  conditions  duplicating  that  of  the 
experiment.  Turbulent  conical  solutions  were  generated  at  x  =  15.2mm  (see 
Figure  3)  and  used  as  initial  data  for  the  PNS  marching  code.  The  present 
computations  used  a  grid  consisting  of  36  circumferential  points  ( A<j)  =  10°) 


8.  Reklis,  R.P. ,  and  Sturek,  W.B. ,  "Surface  Pressure  Measurements  on  Slender 

Bodies  at  Angle  of  Attack  in  Supersonic  Flow",  U.S.  Army  Ballistic 
Research  Laboratory/ ARRADCOM  Memorandum  Report  ARBRL-MR-02876 ,  Aberdeen 
Proving  Ground,  MD,  November  1978.  AD  A064097. 

9.  Kayser,  L.  D. ,  and  Sturek,  W.  B. ,  "Experimental  Measurements  in  the 
Turbulent  Boundary  Layer  of  a  Yawed,  Spinning  Ogive-Cylinder  Body  of 
Revolution  at  Mach  3.0.  Part  II:  Data  Tabulation",  U.S.  Army  Ballistic 
Research  Laboratory /ARRADCOM  Memorandum  Report  ARBRL-MR-02813,  Aberdeen 
Proving  Ground,  MD,  March  1978.  AD  A053458. 

10.  Kayser,  L.  D.,  and  Sturek,  W.  B .,  "Turbulent  Boundary  Layer  Measurements 
on  the  Boattail  Section  of  a  Yawed,  Spinning  Projectile  Shape  at  Mach 
3.0",  U.S.  Army  Ballistic  Research  Laboratory /ARRADCOM  Memorandum  Report 
ARBRL-MR-02880,  Aberdeen  Proving  Ground,  MD,  November  1978.  AD  A065355. 

11.  Sturek,  W.  B.,  "Boundary  Layer  Studies  on  a  Spinning  Cone, "U.S.  Army 
Ballistic  Research  Laboratory /ARRADCOM  Report  BRL-R-1649,  Aberdeen 
Proving  Ground,  MD,  May  1973.  AD  762564. 

12.  Nietubics,  C.  J .,  and  Opalka,  K. ,  "Supersonic  Wind  Tunnel  Measurements  of 
Static  and  Magnus  Aerodynamic  Coefficients  for  Projectile  Shapes  with 
Tangent  and  Secant  Ogive  Noses",  U.S.  Army  Ballistic  Research  Laboratory/ 
ARRADCOM  Memorandum  Report  ARBRL-MR-02991 ,  Aberdeen  Proving  Ground, 

MD,  February  1980.  AD  A083297. 
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and  50  points  radially  between  the  body  and  the  outer  boundary.  Computation 
time  on  the  CDC  7600  computer  is  2.3  sec/step  with  this  size  grid.  A  detailed 
comparison  of  the  PNS  computations  to  these  data  for  a  non-spinning  model  has 
been  accomplished13.  The  comparisons  of  Ref.  13  showed  that  the  PNS  numerical 
technique,  in  which  the  viscous  layer  and  the  inviscid  flow  are  computed 
simultaneously,  yields  a  significantly  improved  agreement  with  experimentally 
measured  wall  pressures  over  the  agreement  obtained  with  inviscid  computations 
for  bodies  with  discontinuities  in  surface  curvature. 

Surface  Pressure.  The  PNS  computations  are  compared  to  experimental 
measurements  an3  to  inviscid  flow  computations  made  using  codes  based  on 
MacCormack's  predictor-corrector  technique  (Figures  4-5).  Longitudinal  sur¬ 
face  pressure  distributions  along  the  windward  and  leeward  rays  are  shown  in 
Figure  4  for  an  angle  of  attack  of  6.3°.  The  PNS  computations  exhibit  better 
agreement  with  experiment  in  the  vicinity  of  the  discontinuities  in  streamwise 
surface  curvature  at  the  ogive-cylinder  and  cyl i nder-boattai 1  junctions  than 
the  inviscid  computations. 

Examples  of  comparisons  of  circumferential  surface  pressure  distributions 
are  shown  for  ot  =  6.3°  in  Figure  5  at  two  longitudinal  stations;  one  on  the 
cylinder  portion  of  the  model  near  the  boattail,  the  second,  midway  on  the 
boattail.  The  comparison  on  the  cylinder  indicates  excellent  agreement 
between  the  PNS  computation  and  experiment  and  the  appearance  of  a  systematic 
discrepancy  between  the  inviscid  computation  and  experiment  for  100°  <  <j>  < 
150°.  This  trend  is  accentuated  for  flow  on  the  boattail. 

Streamwise  Velocity  Profiles.  A  sensitive  test  of  the  accuracy  of  the 
PNS  computational  technique  applied  to  this  flow  is  the  comparison  of  measured 
and  computed  boundary- layer  velocity  profiles.  Such  comparisons  are  shown  in 
Figures  6  and  7  for  two  longitudinal  stations;  station  A  on  the  cylinder  near 
the  boattail,  and  station  B  on  the  boattail  (Figure  3).  Each  figure  shows  the 
velocity  profiles  at  a  particular  longitudinal  station  for  circumferential 
stations  completely  around  the  model  in  30°  increments.  The  nondimensi onal 
streamwise  velocity  component,  u,  is  plotted  versus  physical  distance  y 
measured  radially  from  the  body  surface  in  millimeters,  rather  than  against 
normalized  y/6.  This  method  of  plotting  prevents  scaling  differences  between 
the  computation  and  experiment  from  giving  a  false  comparison. 

Comparisons  for  M  =  3,  a  =  4.2°  and  a  spin  rate  of  333  RPS  are  shown  in 
Figures  6  and  7.  The  agreement  is,  in  general,  excellent.  However,  a  slight 
discrepancy  is  visible  for  the  profiles  just  off  the  lee-side  at  4  =  150°  and 
210°.  This  discrepancy,  which  is  attributed  to  the  formation  of  longitudinal 


13.  Schiff,  L .  B. ,  and  Sturek ,  W.  B .  "Numerical  Simulation  of  Steady 

Supersonic  Floti  Over  an  Ogive-Cylinder-Boattail  Body",  AIAA  Paper  No.  80- 
0066,  18th  Aerospace  Sciences  Meeting,  January  1980. 
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vortices,  becomes  more  pronounced  at  increased  angles  of  attack  and  is  discus¬ 
sed  in  more  detail  in  Ref.  13.  For  the  purpose  of  this  study  -  which  is  to 
determine  Magnus  effects  at  small  angle  of  attack  -  the  agreement  between 

computational  and  experimental  velocity  profiles  is  considered  to  be  very 
satisfactory. 

Magnus  Force.  The  Magnus  effect  is  strongly  dependent  on  the  circumfer¬ 
ential  pressure  distributions  and  the  asymmetries  about  the  pitch  plane  that 
occur  when  the  model  is  spinning.  Thus  the  PNS  technique  is  particularly 

attractive  for  computations  of  the  Magnus  effect  since,  as  seen  in  Figures  4- 
5,  the  accuracy  for  the  circumferential  and  longitudinal  distributions  of  wall 
pressure  are  much  improved  over  that  obtained  using  inviscid  techniques.  This 
improvement  was  most  significant  for  the  flow  over  the  boattail.  In  order  to 
test  the  accuracy  of  the  PNS  code  for  computing  Magnus  effects,  several  test 
cases  have  been  run  for  flow  conditions  where  experimental  measurements  of 

Magnus  forces  and  moments  have  been  accomplished11’ 12  .  See  Figure  8  for  the 
sign  convention  used  in  evaluating  the  aerodynamic  forces  and  moments. 

The  results  of  these  computations  for  the  10°  cone  model  are  shown  in 
Figure  9  and  for  the  ogive-cylinder  and  ogi ve-cyl i nder-boattai 1  models  in 
Figures  10  and  11,  respectively.  The  PNS  computations  are  compared  to 

i nvi sci d- boundary  layer  results  and  to  the  experimental  measurements  as 
indicated.  Tabulated  results  are  presented  in  Tables  1,  2  and  3. 

In  Figures  9,  10  and  11,  the  individual  components  of  the  Magnus  effect 
are  plotted  versus  axial  position  for  Mach  =  3,  a  =  2°,  and  ft  =  333  RPS. 
Comparisons  are  shown  of  PNS  and  INV-BL  results  for  the  individual  components 
of  the  Magnus  effect.  Excellent  agreement  for  CT(j)  (circumferential  wall  shear 

component)  is  indicated.  Fair  agreement  is  indicated  for  Cpw  (wall  pressure 

component  of  Magnus)  between  the  two  computational  techniques.  For  the  INV-BL 
technique,  the  component  plotted  as  Cpw  is  the  sum  of  the  wall  pressure  and 

centrifugal  pressure  gradient  contributions.  This  quantity  is  equivalent  to 
the  surface  pressure  contribution  determined  by  the  PNS  code  since  the  viscous 
and  inviscid  flow  are  computed  simultaneously. 

The  lines  in  Figures  9-11  identified  as  Cy  represent  the  total  Magnus 

effect.  For  the  cone,  both  the  PNS  and  INV-BL  techniques  achieve  good  agree¬ 
ment  with  the  experiment.  Considering  the  small  magnitude  of  the  Magnus 
force,  the  agreement  with  experiment  shown  in  Figure  10  for  the  ogive-cylinder 
model  is  considered  to  be  quite  good  for  the  INV-BL  technique;  however  the  PNS 
result  is  obviously  better.  This  result  is  also  considered  to  be  verification 
of  the  INV-BL  concept  at  small  angle  of  attack.  The  results  shown  in  Figure 
11  show  a  dramatic  superiority  of  the  PNS  technique  compared  to  the  INV-BL 
technique  for  the  ogive-cylinder  boattail  model  where  the  INV-BL  technique 
greatly  exaggerates  the  effect  of  the  flow  over  the  boattail. 

Comparisons  of  the  PNS  computations  to  experiment  for  additional  Mach 
number  flow  conditions  (2  <  M  <  4)  for  the  ogive-cylinder  and  ogi ve-cyl i nder- 
boattai  1  models  are  shown  for  Maqnus  force  coefficient  in  Figure  12  and  for 
the  slope  of  the  Magnus  moment  coefficient  in  Figure  13.  Excellent  agreement 
is  indicated  for  the  Magnus  force  coefficient  in  Figure  12  for  magnitude  and 
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or  trend  with  Mach  number  and  body  configuration.  The  agreement  shown  for 
the  slope  of  the  Magnus  moment  coefficient  in  Figure  13  is  less  satisfactory; 
however,  the  trend  with  Mach  number  and  body  configuration  is  accurately 
predicted.  The  absolute  value  of  the  prediction  is  well  within  that  required 
for  parametric  design  studies. 

Computations  have  also  been  performed  using  the  PNS  code  for  laminar 
boundary  layer  flow  over  a  10°  cone  at  M  =  4,  a  =  2° ,  S2  =  500  RPS, 

r”/To  =  °*2T  A  constant  wall  temperature  boundary  condition  was  included  in 
the  PNS  code  for  this  computation.  This  is  the  case  for  which  a  reversed 
Magnus  effect  was  reported  in  Ref.  5.  The  results  of  these  computations  are 
shown  in  Figures  14  and  15.  Individual  components  of  the  side  force  coef¬ 
ficient  due  to  the  Magnus  effect  obtained  with  both  computational  methods  are 
shown  in  Figure  14.  The  results  differ  considerably.  In  contrast  to  the 
r  esults  of  Ref.  5,  which  show  a  reversal  of  the  wall  pressure  component  and 
total  Magnus  force  with  increasing  body  length,  the  present  PNS  results  vary 
monotonical ly  with  body  length.  (Note  that  in  Ref.  5  values  of  CTX  and  Cy  are 

given  only  at  x  =  304.8  mm  (1.0  ft).  The  intermediate  values  shown  in  Figure 
14  were  obtained  by  assuming  a  linear  variation  of  CTX  with  body  length  and 

summing  the  individual  components  CTX,  CT(j>,  and  Cpw  to  obtain  Cy.)  Additional 

computations  were  made  for  a  wide  range  of  grid  configurations  and  the  results 
evidenced  no  reversal  of  the  Magnus  force.  The  monotonic  behavior  of  the 
present  PNS  results  follows  the  classical  trend,  and  is  similar  to  that  shown 
above  for  the  turbulent  case. 

The  results  of  an  investigation  of  the  effect  of  wall  temperature  is 
shown  in  Figure  15,  where  computations  for  a  cold  wall  and  adiabatic  wall 
boundary  conditions  are  compared.  No  tendency  for  a  reversal  of  the  Magnus 
force  was  obtained.  The  effect  of  the  cold  wall  is  to  reduce  the  magnitude  of 
the  Magnus  force. 


C.  Pitch  Plane  Aerodynamics 

Examples  of  computed  values  for  the  slope  of  the  pitching  moment  coef¬ 
ficient,  Cm  ,  and  center  of  pressure,  CP,  are  shown  in  Figures  16  and  17 

a 

respecti vely,  compared  to  experiment  for  the  ogi ve-cyl inder  and  ogive- 
cylinder-boattail  shapes.  The  agreement  with  experiment  is  excellent  for  the 
full  range  of  Mach  number,  2  <  M  <  4. 


IV.  CONCLUDING  REMARKS 

A  numerical  computational  study  has  been  described  in  which  the  PNS 
marching  code  recently  developed  by  Schiff  and  Steger  has  been  used  to  compute 
Magnus  effects  for  spinning,  slender  bodies.  The  PNS  computations  have  been 
compared  to  INV-BL  computations  and  to  experimental  measurements.  The  compar- 
isons  have  shown  that  the  PNS  code  yields  excellent  agreement  with  the  experi¬ 
mental  data.  These  results  indicate  that  the  PNS  numerical  technique,  which 
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computes  the  inviscid  flow  and  the  viscous  layer  simultaneously,  closely 
models  the  physics  of  the  actual  flow  over  the  body  including  discontinuities 
in  surface  curvature.  This  feature  has  been  shown  to  be  of  particular  signif¬ 
icance  for  computations  of  the  flow  over  a  boattailed  afterbody. 

The  results  also  indicate  verification  of  the  suitability  of  the  INV-BL 
technique  for  cone  and  ogive-cylinder  bodies  at  small  angle  of  attack  in 
contradiction  to  a  previous  result  for  laminar  viscous  flow. 

The  comparisons  of  this  paper  represent  the  first  comprehensive  tests  of 
PNS  computations  for  turbulent,  viscous  flow  over  a  spinning  ogi ve-cyl i nder- 
boattail  model  and  indicate  that  the  PNS  code  reported  here  represents  a 
viable  computational  tool  for  predicting  Magnus  effects  for  spinning  pro¬ 
jectiles  at  small  angle  of  attack. 
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PNS 

INV-BL 

EXP 


Table  1.  Magnus  Force  Components  for  10°  Cone- 
■  M  =3:T  a  =■  2°,  a  =  333  RPS 

^TX  ^T<j>  ^  Ap  ■$..  Cpw  Cy 

-.769x10” 6  .772x 10" **  -.933xl0”7  -.244xl0"3  -.168xl0”3 

•636xl0"6  .792xl0”4  -.170xl0'3  -.108xl0"3  -.198xl0"3 

-.170x10" 3 


PNS 

INV-BL 

EXP 


Table  2.  Magnus  Force  Components  for  Ogive-Cylinder  Body- 
,M  -3,  a  =2°,  & =  333  RPS 

CTX  ^T(j)  CAp  ^pw  CY 


-.102xl0”5  .215xl0”3  -.842xl0"6  -.253xl0"2  -.00232 

.570x10  6  .197xl0”3  -.161xl0"2  -.242xl0"2  -.00383 

-.00250 


Table  3.  Magnus  Force  Components  for  Ogive-Cylinder-Boattail  Body 

M  =  3,  a  =  2°,  fl  =  333  RPS 


CTX 

CT(J) 

CAp 

^pw 

Cy 

PNS 

-.147xl0”5 

•180xl0”3 

-.801xl0”6 

-.278xl0”2 

-.00259 

INV-BL 

•225xl0”5 

•175xl0”3 

-.190xl0”2 

-  .427xl0”2 

-.00599 

EXP 

-.00300 
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Figure  1.  Coordinates  and  Notation 


*. 


20 


u  =  a 
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Figure  2*  Subsonic  Layer  Approximation 
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fl,B  BOUNDARY  LAYER  SURVEY  STATIONS 
C  CONICAL  STARTING  SOLUTION  STATION 

Figure  3.  Model  Configuration  for  Detailed  Flow  Field  Studies 
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Figure  4.  Axial  Surface  Pressure  Distribution  for  Ogive-Cyl inder-Boattail 
Body;  M  =  3,  gt  =  6.3°,  fi  =  0  RPS,  Re  =  2.13  x  107/m 
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_  pNNVsISCl0} COMPUTATIONS 

o  EXPERIMENT  REF  8 


Figure  5.  Circumferential  Surface  Pressure  Distributions  for 
Ogive-Cyl inder-Boattail  Body;  M  =  3,  a  =  6.3°, 
a  =  0  RPS,  Re  =  2.13  x  io7/m 
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Figure  6.  Boundary-Layer  Velocity  Profiles  for  Ogive-Cyl inder-Boattail 
Body;  M=3,a=4.2°,n=  333  RPS,  Re  =  2.13  x  107/m, 
x  =  254mm,  Cylinder 
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PNS  COMP 
EXP,  REF  10 


Figure  7.  Boundary-Layer  Velocity  Profiles  for  Ogive-Cyl inder-Boattail 
Body;  M  =  3,  a  =  4.2°,  n  =  333  RPS,  Re  =  2.13  *  107/m, 
x  =  324mm,  Boattail 
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NORMAL  FORCE 
CN 


MAGNUS  (SIDE)  FORCE 


Figure  8.  Sign  Convention  for  Aerodynamic  Forces 
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Figure  9. 


Magnus  Force  Coefficient  for  10° 
Q  =  333  RPS,  Re  =  2.352  x  107/m, 


Cone;  M  =  3,  a  =  2°, 
Turbulent  Boundary-Layer 


28 


*  10*  SIDE  FORCE  COEFFICIENT 


Figure  10.  Magnus  Force  Coefficient  for  Ogive-Cylinder  Body;  M  =  3,  a  =  2° 
-  333  RPS,  Re  -  2.11  x  10^/m,  Turbulent  Boundary-Layer 
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103,  SIDE  FORCE  COEFFICIENT 


Figure  11.  Magnus  Force  Coefficient  for  Ogi ve-Cyl inder-Boattail  Body; 

M  =  3,  a  =  2°,  ft  =  333  RPS,  Re  =  2.11  x  107/m,  Turbulent 
Boundary-Layer 
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Figure  12b.  Magnus  Force  Coefficient  for  Ogive-Cyl inder-Boattail 
Body;  a  =  2°,  n  =  333  RPS,  Turbulent  Boundary-Layer 
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Figure  13a.  Magnus  Moment  Coefficient  for  Ogive-Cylinder 
Body;  Turbulent  Boundary  Layer 


Figure  13b.  Magnus  Moment  Coefficient  for  Ogi ve-Cyl inder-Boattail 
Body;  Turbulent  Boundary-Layer 
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Figure  14.  Magnus  Force  Coefficient  for  10°  Cone;  M  =  4,  a  =  2°, 
ft  =  500  RPS,  Tw/To  =  0.24,  Re  =  9.22  x  106/m,  Laminar 
Boundary-Layer 
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lO3,  SIDE  FORCE  COEFFICIENT 
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Figure  15. 

Magnus  Force  Coefficient  for  10°  Cone;  M  =  4,  a  =  2°, 
n  =  500  RPS,  Re  =  9.22  *  106/m,  Laminar  Boundary-Layer 

34 


Figuer  16a.  Pitching  Moment  Coefficient  for  Ogive-Cylinder 
Body;  Turbulent  Boundary-Layer 


Figure  16b.  Pitching  Moment  Coefficient  for  Ogi ve-Cyl inder- 
Boattail  Body;  Turbulent  Boundary-Layer 
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Figure  17a.  Center  of  Pressure  for  Ogive-Cylinder 
Body;  Turbulent  Boundary-Layer 


Figure  17b.  Center  of  Pressure  for  Ogi ve-Cyl inder- 
Boattail  Body;  Turbulent  Boundary-Layer 
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LIST  OF  SYMBOLS 


a  speed  of  sound 

Cp  specific  heat  at  constant  pressure 

CP  center  of  pressure  for  normal  force 

Cm  pitching  moment  coefficient 

Cm  dCm/da  »  slope  of  pitching  moment  coefficient  evaluated  at  a  =  0 

a 

Cn  Magnus  (yawing)  moment  coefficient 

C1  d2Cn/[d(-^)  *da]  ,  slope  of  Magnus  moment  coefficient  evaluated  at 

pa  0  =  0,  a  =  0 


CY 

C 


pw 


'Ap 


xx 


'x<i> 


D 

e 

ei 

A  AAA 

,F,G»q 


L 

M 

P 

Pr 

Re 

A 

Re 


Magnus  (side)  force  coefficient 

wall  surface  pressure  component  of  Maqnus  force  coefficient 

centrifugal  pressure  gradient  component  of  Magnus  force  coefficient 

longitudinal  velocity  wall  shear  component  of  Magnus  force 
coeffi cient 

circumferential  velocity  wall  shear  component  of  Magnus  force 
coefficient 

diameter  of  model 

total  energy  per  unit  volume  of  fluid,  normalized  by  p^a^2 
internal  energy,  normalized  by  a^2 

flux  vectors  of  transformed  gasdynamic  equation  (Eq.  2) 

Jacobian  of  transformation  between  physical  and  computational 
coordinates 

reference  length 

Mach  number 

pressure,  normalized  by  p^a,,2 
Prandtl  number,  y  c  /< 

7  roo  p'  00 

Reynolds  number,  p  U  L/y 
00  00  '  00 

Reynolds  number  (Eq.  1),  P^a^L/y^ 
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s 


u,v,w 


U,V,W 

x,y,z 

a 

Y 

K 


P 


C,n,c 


p 

<p 

S2 


LIST  OF  SYMBOLS 
(Continued) 

viscous  flux  vector  (Eq.  6) 

Cartesian  velocity  components  along  the  x,  y,  z  axis,  respectively, 
normalized  by  a 

contravariant  velocity  components  (Eq.  3) 
physical  Cartesian  coordinate  axes  (Fig.  1) 
angle  of  attack 
ratio  of  specific  heats 

coefficient  of  thermal  conductivity,  normalized  by  free- stream 
value  < 

oo 

coefficient  of  viscosity,  normalized  by  free-stream  value  u 

oo 

computational  coordinates  in  the  axial,  circumferential,  and  radial 
directions  (Fig.  1) 

density,  normalized  by  free-stream  density 
circumferential  angle  (Figure  1) 
spin  rate  about  body  axis,  RPS 


Subscri pts 

“  free-stream  conditions 

w  body  surface  values 

x  based  on  axial  distance  from  nose 
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